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Abstract. We present a study of radial density structure of the star forming globule, DCld303.8-14.2 (DC303), and a non-star 
forming globule, Thumbprint Nebula (TPN), using near-infrared data taken with the ISAAC instrument on the Very Large 
Telescope. We derive the extinction through the globules using the color excess technique and examine the radial density 
distribution using Bonnor-Ebert and power-law models. The two globules have significantly different density structures. The 
extinction profile of DC303 is best fitted with a single power-law with an exponent p = 2.29 ± 0.08. An unstable Bonnor-Ebert 
model with a dimensionless parameter £ max = 23 ± 3 provides equally good fit to data. The extinction profile of TPN flattens at 
small radii, making the profile significantly different from the profile of DC303. We are unable to fit the Bonnor-Ebert model 
for TPN in a robust manner, but derive the lower limit £ max > 8 for the dimensionless outer edge. The density profile derived 
for TPN is typical compared to recently observed pre-protostellar globules, with high £ max value which could be interpreted as 
the presence of significant additional support or very slow contraction. 

Key words. Stars: formation - ISM: clouds - ISM: individual objects: DCW303.8-14.2 - ISM individual objects: Thumbprint 
Nebula 



1. Introduction 

Small-sized clouds of cold interstellar matter, the globules, are 
often referred to as ideal places to study processes and phenom- 
ena connected to low-mass star formation. The ideality arises 
from the apparent simplicity of these objects: they are nearly 
spherical, isolated, slowly evolving clouds and only a small 
number of stars, if any, can be formed in their interiors. The 
sizes of typical globules are in the order of tenths of parsecs, 
their masses range from ~ 1 M ... 20 M , and the tempera- 
ture in them have been observed to be around 10 K. Roughly 
one third of the globules contain point sources, indicating the 
presence of a newly born star (Lee & Myers 1999). 

Send offprint requests to: J. Kainulainen e-mail: 
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Once a globule has been formed in a molecular cloud it may 
evolve towards a point where the gravitational force inevitably 
takes over and collapse ensues. In theories of star formation the 
properties of the entire collapse process depend on the char- 
acteristics of the cloud at the onset of the collapse, i.e. at the 
initial stage of the collapse. Factors such as mass distribution, 
temperature, velocity and magnetic fields at initial stage pre- 
dict the future evolution of the collapsing core. Some of the 
factors, like temperature and large scale velocity fields, can be 
well traced with observations of molecular lines. Further, the 
significance of some factors, e.g. magnetic fields, are not so 
well understood due to observational difficulties. 

In particular, theoretical star formation scenarios predict the 
density distribution for the core at the initial stage of the col- 
lapse. The form of the density distribution is linked to the rate 
at which a protostar accretes mass from the surrounding en- 
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velope and further to the timescale of the collapse. It is thus an 
essential factor considering the process of star formation. In the 
classical 'inside-out' collapse scenario of Shu (Shu 1 19771 S77 
hereafter) the initial stage of core collapse equals to the singular 
isothermal sphere with a steep density distribution p(r) ~ r~ 2 . 
After this configuration is established the collapse begins in 
the center of the core expanding "inside-out" as a shock-front. 
The mass accretion rate in S77 is constant throughout the col- 
lapse for all radii. Some other models, such as those of Foster 
& Chevalier ( 1993 ) who take a marginally stable Bonnor-Ebert 
sphere as an initial condition, or Ciolek & Mouschovias ( 1994) 
and Basu & Mouschovias (119941) who take ambipolar diffusion 
and magnetic fields into account, assume flat inner density pro- 
files for the initial stage of the collapse. In those models the 
accretion rate is a function of both radius and time. 

The observations of pre-stellar cores having both p(r) ~ 
r~ 2 and flat inner profiles indeed have been made using sub- 
millimeter dust continuum observations (e.g. Ward-Thompson 
et al. [19941 [19991 Chandler & R icher |2000l Shirley I2U001 
Johnstone et al. 120001 I200T1 [20051 Evans et al. I200TI Kirk et 
al. 120051 120061 ) . For some individual globules the mass dis- 
tribution have been obtained via molecular line observations 
and chemical modeling (e.g. Zhou et al. 119931 van der Tak et 
al. |2005 ). However, derivation of the density distribution us- 
ing these methods, i.e. interpreting the observations in terms of 
density, suffers from the unknown properties of the selected 
model: excitation conditions, possibly large optical depths, 
time dependent chemistry or unknown temperature and emis- 
sivity of the dust grains. 

The extinction of background starlight caused by dust par- 
ticles provides a more direct method than those above for de- 
termination of the density distribution. The colors of the stars 
shining through the cloud are measured at two or more bands 
and the colors are compared to the colors of stars in a nearby 
reference field. The near-infrared is particularly suitable for the 
observations as the extinction is about one tenth of the extinc- 
tion in visual wavelengths and cores therefore more transpar- 
ent. This method has been successfully used in the context 
of dark globules, leading for example to the suggestion that 
pre-protostellar cores could be well described by pressure con- 
fined Bonnor-Ebert models close to the hydrostatic equilibrium 
(Alves et al. [2001). Moreover, protostellar cores which show 
strong evidence of star formation activity have been observed 
to have density distributions corresponding to highly unstable 
Bonnor-Ebert spheres or steep power-laws with exponent p > 2 
(Harvey et al. l^TXTTl Harvey et al. 120031) . 

The observations of density distributions from both pre- 
protostellar and protostellar cores can be used to seek the char- 
acteristics of different star formation scenarios. In this paper we 
exploit the near-infrared color excess technique to study the ra- 
dial density distributions of two cloud cores: we determine the 
mass distribution of a core that has very recently formed a star, 
and a core which is regarded to be quiescent, close to hydro- 
static equilibrium and thus possibly close to the initial stage of 
star formation. For this purpose we have selected two globules, 
namely DCW303. 8-14.2 (or Sandquist 160) and Thumbprint 
Nebula, and conducted deep near-infrared JHK surveys to- 
wards them. 



The two globules are introduced in §J2] and in ^3] we de- 
scribe the observations, data reduction and the photometry. The 
radial density distributions of the globules are presented and 
discussed in ^4] and |J5] The conclusions of the work are pre- 
sented in $6] 

2. DCId303.8-14.2 and Thumbprint nebula 

The globule DCW303.8-14.2 (DC303) and Thumbprint Nebula 
(TPN) are both members of the nearby, quite extensively stud- 
ied molecular cloud complex of Chamaeleon (Cha). They are 
located close to the clouds Cha II and Cha III, from which Cha 
II is an active star forming cloud with tens of Young Stellar 
Object (YSO) candidates within its boundaries and Cha III is 
a quiescent cloud. The two globules, DC303 and TPN, were 
initially selected for this study due to their similar morpholog- 
ical appearance and accessible location in a well-studied cloud 
complex. 

Fig. [T] shows the IRAS 100 //m map of Cha II/III region 
and the Digitized Sky Survey blue plate images of both glob- 
ules. The globules are separated by about 2°. Both cores are 
clearly visible in the 100 /vm map and DC303 is actually one 
of the brightest 100 /mi sources in Cha II after the central star 
forming region. The globules show very similar morphology in 
the DSS blue plate images (Fig. []]): the surface brightness of 
scattered light forms a bright rim around the globule with max- 
imum intensity on the northern side, which is the direction of 
the galactic plane. Inside the rim the surface brightness drops 
forming a dark center for the globules. 

DC303 harbors an IRAS point source IRAS 13036- 
7644. The millimeter wavelength emission lines of 12 CO(/=l- 

0) , CS(7=2-1 ) and (7=5-4), HCO + (7=1-0), HCN(/=l-0) 
by Lehtinen (1997) and emission lines of H2CO(2i2-ln), 
C 3 H 2 (3i2-2 2 i) and CH 3 OH(3 -2 A + ) by Mardones et al. 
(1997J) have revealed a low- velocity outflow driven by the 
IRAS source. Also, some molecules show blue-skewed line 
profiles indicating the possibility of collapse motions. On the 
other hand, some lines such as HCO + show red-skewed pro- 
files, and radiative transfer modeling would be required to 
further interpret the molecular line observations. The spectral 
energy distribution of IRAS 13036-7644 is typical of a very 
young protostar. Based on the characteristics of the SED and 
qualitative arguments about the molecular line profiles, the 
source is placed in the midway of Classes and I in the tra- 
ditional classification scheme of YSOs (Lehtinen et al. [2003b, 
120531 see e.g. Lada & Wilking [19841 Adams et al. [19871 
Wilking et al. 1989 for the classification scheme). 

The Thumbprint Nebula is an apparently quiescent glob- 
ule with no detected infall/outflow motions or embedded point 
sources. It has been studied previously in near-infrared to con- 
strain dust model parameters (Lehtinen et al. 119961) and in far- 
infrared to derive the dust temperature (Lehtinen et al. 1995; 
Lehtinen et al. 1 1 998b . The millimeter wavelength molecular 
lines of 12 CO(/=1-0) and (7=2-1), 13 CO(/=1-0) and (7=2- 

1) , C 17 O(/=l-0), CS(7=2-1) and HCN(/=l-0) have been ob- 
served towards the globule, from which a very low kinetic tem- 
perature Tkin ~ 6.6 K was derived. The nebula also seems to be 
close to the thermodynamic equilibrium with E po ,IEkm ~ 0.9 
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(Lehtinen et al. 1 1 9951 >. Lehtinen et al. also find that the energy 
balance of the globule is explained solely by the external radi- 
ation from the interstellar radiation field. 

3. Observations and data reduction 

The observations were made in service mode with the ISAAC 
(Infrared Spectrometer And Array Camera, Moorwood et al. 
1 19981 1 at VLT (Very Large Telescope) during January-April 
2002, in J, H and K s bands. 

In addition to the measurement of extinction, the objective 
of the observations was to detect the weak surface brightness 
of the scattered light. Because the extent of the scattered light 
was expected to be about the same as the field of view of a sin- 
gle image frame, the observations were made in an OFF-ON- 
OFF. . .-ON-OFF sequence, where the OFF-position is outside 
the cloud. Individual ON and OFF frames were taken as 30 
sec or 1 min exposures. The OFF-positions were then used to 
derive the sky background image. We observed several ON po- 
sitions for both clouds to fully map the clouds and to have some 
background sky around the clouds. The positions of the ON and 
OFF fields are shown in Fig.Q]for both globules. 

The basic data reduction, such as sky substraction and flat- 
fielding, has been performed with the xdimsum package within 
IRAF. To correct for the variability of the sky background, 
the 3-6 OFF-positions closest (in time) to each dithered ON- 
position were combined and substracted. These sky-substracted 
ON-frames were then shifted and added to create a total ON- 
position image of typically 20-30 min total exposure time. 

Then the SExtractor (Bertin & Arnouts 1996) package was 
used to find point sources in the images. The magnitudes of the 
stars have been derived with the IDL version of the DAOPHOT 
aperture photometry package (http://idlastro.gsfc.nasa.gov). 
The limitting magnitudes, defined as magnitudes whose pho- 
tometric error is equal to 0.1 mag, are: J - 20.8, H - 19.2 
and K s = 18.8 mags for DC303 and J = 20.7, H = 19.0 and 
Ks = 18.7 for TPN. Typical errors of J, H and Ks band magni- 
tudes are 0.05 mag, 0.08 mag and 0.08 mag, respectively. 

The photometric calibration was derived using stars from 
the Persson faint NIR standard list (Persson et al. 1998). 
Several photometric nights during the observations were used 
to fix the photometric scale, and the zero-points for the non- 
photometric nights were derived by utilizing overlapping fields 
from different nights. Overall photometry is uniform over the 
whole field within 5%, and the JHKs magnitudes were also 
checked to be consistent with 2MASS photometry within few 
percents. The photometry of all detected sources will be made 
available online at the CDS. 

4. Results 

4.1. NIR images of DC303 

The reduced center frames of DC303 observed with ISAAC at 
Ks, H and J bands are shown in Fig. [2] The images show two 
distinctive features: the bright circle of scattered light forming a 
circular halo around the center of the core and a number of non- 



pointlike sources close to the location of the embedded IRAS 
source (see also Lehtinen et al. 2003a). 

The scattered light has been commonly observed towards 
dark clouds at optical wavelenghts, but applications to near- 
infrared wavelengths have been few. Dedicated observations 
of scattered light in near-infrared have been made to con- 
strain dust models in dark clouds (Lehtinen & Mattila 1996; 
Nakajima et al. 120031 ). but the phenomenon has not been ex- 
ploited in the context of tracing the mass distributions. Several 
studies in which NIR imaging of dark cores have been con- 
ducted have, however, reported the presence of diffuse struc- 
tures (Lehtinen et al. I2003al Kandori et al. 120051 ). Recently, 
Foster & Goodman (2006 ) and Padoan et al. (120061 ) studied the 
feasibility of using measurements of scattered light in near in- 
frared to derive the column density distribution of a dark cloud 
and concluded that it can be used, when complemented with 
a proper radiative transfer modeling, as a new high-resolution 
tool. The measurements of scattered light seem to be particu- 
larly well suited to determine the density distribution of more 
diffuse clouds, in which the maximum extinction remains be- 
low ~ 20 mag (Juvela et al. 2006). Although the extinction in 
the center of DC303 is clearly higher, it offers an interesting 
ground to implement the treatment into the case of dark glob- 
ules. Our aim is to model the scattered light to constrain the 
mass distribution and study the dust properties in DC303 in a 
forthcoming paper. The extended structures close to the IRAS 
source resemble the appearance of Herbig-Haro objects often 
seen associated with outflows from young stars. We have per- 
formed near-infrared long-slit spectroscopy of these objects to 
verify their nature. 

4.2. NIR colors and extinction 

We use the H and K magnitudes of the observed stars to derive 
the near-infrared extinction in the globules. In particular, we 
use the color excess technique described by Lada et al. ( 1994, 
"NICE"). The color excess along the line of sight towards a star 
is related to the observed magnitudes through the equation: 

E(H - K) = (H - K) obs - (H - K) int , (1) 

where (H - K) mi is the intrinsic color of the star. In the NICE 
method, this equation is solved by assuming that the intrinsic 
colors of background stars are, on the average, similar to the 
colors of stars observed on a nearby field which is free from 
extinction. Furthermore, the color excesses can be converted to 
visual extinction by assuming a certain form of the reddening 
law, i.e. the ratio between color excess and visual extinction. 
Using e.g. the often quoted reddening law of Rieke & Lebofsky 
(11985b we obtain: 

A y = 15.87 x E(H-K) (2) 

This conversion is not necessary for the derivation of radial 
density distributions in Sec. 14.31 but is provided to have a com- 
mon ground for comparison with the results of other papers. 

The average intrinsic colors calculated from the reference 
fields (Fig. [Q) are: (H - K s ) = 0.14 + 0.13 for DC303 and 
(H - K s )o = 0.13 + 0.09 for TPN. Applying Eqs. [T] and |2 
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Fig. 1. Left: IRAS 100 fim emission map of the Chamaeleon II-III region. The white frames show the positions of the DSS images 
on the right. Center: The Digital Sky Survey blue band image of the DC303 region. The frames observed with the VLT/ISAAC 
are marked with white rectangles. The position of the point source IRAS 13036-7644 is marked with an asterisk. Right: The 
same for the Thumbprint Nebula. 




Fig. 2. The K s , H and J band images of DC303. The asterisk marks the position of the IRAS 13036-7644. 



and using these intrinsic colors, an estimate is obtained for ex- 
tinction towards each observed star in the globule fields. We 
assume that all the stars are background stars (i.e. there are no 
foreground stars). To produce a spatial distribution from the 
discrete data, the color excess values are then convolved using 
a gaussian function with 40" and 30" FWHM for DC303 and 
TPN, respectively. Each pixel in the final map is required to 
have at least one star within the FWHM radius from the pixel 
center to make the extinction calculation applicable. If there 
are no stars within that range, the extinction in the pixel is con- 
sidered too high to be determined using this method. The ex- 
tinction maps resulting from this treatment are shown in Fig. 

m 

The extinction map of DC303 reveals a high extinction gra- 
dient towards the position of the IRAS source, forming a sym- 
metric core projected on the plane of the sky. There is a region 
of about 20" in radius in which no H band stars are detected. 
Since at least one star within the FWHM range is required, this 
results in a few pixels in which the extinction is too high to 
be determined. These pixels are marked white on the extinc- 
tion map. The maximum extinction in the map, Ay = 28 m , is 
located at the pixel closest to the IRAS source. 



The extinction in TPN is significantly lower than in DC303 
(Fig. [3] right panel). The maximum, Ay = 10.2™, is located 
at the coordinates (ra, dec)=(12 ;, h44 m 55 s .6, -78 d 48 m W). The 
core appears to be slightly elongated in northwest-southeast 
direction and the gradient of extinction on the northwest side 
of the maximum is clearly higher than on the southeast side. 
The size of the core is quite similar with DC303, both having 
Ay > 3'" extending approximately 2' in radius. 

4.3. Models of radial density distribution 

Fig. |4] shows the (H - Ks) colors of all the observed stars as 
a function of distance from the center of the cloud. In DC303 
the position of the IRAS source has been selected to represent 
the center point of the core. In TPN the center point is chosen 
to be the point of maximum extinction of the Ay map. The 
solid lines in the figures are the gaussian convolved averages 
of the colors in annuli 20" and 15" wide, for DC303 and TPN 
respectively. The dashed lines give the standard deviation of the 
colors. In DC303, the star closest to the IRAS source lies at the 
angular distance of r ~ 25" and possesses the highest (H - Ks) 
color in the sample, i.e. (H - Ks) = 3.08 corresponding to 
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Fig. 3. Left: The visual extinction map of DC303 derived from the VLT/IS AAC observations using the NICE color excess 
method. The FWHM resolution of the map is 40". The white pixels result from not having a single star within the FWHM range 
from the pixel center. The plus sign marks the position of IRAS 13036-7644. Contours start from 3.5 mags and the the step is 2.5 
mags. Right: The same for TPN. The FWHM resolution of the map is 30". Contours start from 1 mag and the the step is 1 mag. 




Fig. 4. Left: The H - K color as a function of the distance from the IRAS point source. The solid line is gaussian convolved 
average of the colors and dashed lines the standard deviation. Right: The same for TPN. 
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about A v ~ 47'". In TPN the star closest to the center point is 
located at the distance of 5" and has the color (H - K§) = 0.93 
corresponding to the visual extinction of Ay as 13™. 

In the following we model the observed column density 
profile using two models for radial density distribution: the 
spherically symmetric inside-out collapse scenario of S77, 
which essentially means the power-law density distribution, 
and the Bonnor-Ebert model (Ebert [19551 Bonnor [19561) . The 
radial density profiles provided by the models are converted 
to column density profiles and compared to the normalized 
(H - K&) data values (Fig. [4j. To evaluate the goodness of fit 
the chi-square value between the model profile and the data is 
used: 



2 = — y f 

dof L 



X = 



E(H - K s )f s - E(H - K s y! 

~.obs 



model 



(3) 



The dof stands for the degrees of freedom in the fit, and cr° bs is 
the uncertainty related to each observation. The uncertainty is 
evaluated as a quadratic sum of the photometric uncertainty of 
the observation, typically 0.08 mag, and the uncertainty in the 
intrinsic color: 



+ o- bg . 



(4) 



4.3.1. Inside-out collapse 



The basic scenario of spherical symmetric collapse of cloud 
cores, the "inside-out" collapse, is described by S77. The 
inside-out collapse begins from the initial stage where the 
core is relaxed to a balance between gravitation and pressure, 
and the density distribution of a singular isothermal sphere, 
p(r) ~ r~ 2 , is established. The collapse begins in the center 
of the cloud, and propagates then outwards as a shock front 
with a speed equivalent to the speed of sound. While the den- 
sity profile in the outer envelope remains at ~ r~ 2 , the density 
profile of the collapsing region approaches p ~ r~ 3 ^ 2 , indicative 
of free-fall motion. 

The normalization of the density distribution at the static 
initial condition is determined by the effective sound of speed, 
a = (fcr/yuw//) 1 ^ 2 , where T is the kinetic temperature and p the 
mean molecular weight. The complete equation for the density 
distribution has the form: 



Pstatici 1 ") 



(5) 



2nGr 2 ' 

The core of DC303 appears to be very symmetric to the ex- 
tent of about 180". The globule itself is extended mainly west- 
wards from this densest core. Only the stars closer than 180" 
are included in the fit to avoid the contamination from this more 
extendend part. 

If the observed extinction profile of DC303 is assumed to 
originate merely from the core itself, the best-fit power-law has 
the exponent p = 2.03 +0.04 (x 2 = 1 .58). If an additional "fore- 
screen" component is included as a free parameter to mimic a 
more diffuse layer surrounding the core, the fit quality improves 
resulting in a slightly higher exponent p = 2.29 + 0.08. In the 
best fit the forescreen component has the value E(H-K) = 0.09 
(X 1 - 1.54). Considering the extinction extending westwards 



from the core the inclusion of the additional component seems 
well justified. The power-law fit with a forescreen component 
is illustrated in Fig. [5] The fit with a forescreen extinction is not 
significantly dependent of the selected radius of the core. The 
profile in Fig.[5]is calculated with R 01 „ = 250". 

We find that the scaling set by the static condition of S77 
(Eq. |5]l matches the observed scale of the extinction well. 
We have used the temperature T = 1 1 K evaluated from the 
C ls O data (Kainulainen 2004), the reddening law of Rieke & 
Lebofsky ( 1985) and standard gas-to-dust ratio of Bohlin et al. 
(H9781 see Eq.|9). 

The extinction profile of TPN at r > 10" is well fitted with 
quite a flat power-law having exponent p = 1.19 + 0.09 (y 2 = 
1.13). However, the three datapoints with r < 10" show clearly 
a turn-over in the profile, and thus a power-law-like model is 
not really reasonable for TPN. 

4.3.2. Bonnor-Ebert models 

Recently, pressure-confined, isothermal gas spheres known as 
"Bonnor-Ebert spheres" (BE) have been succesfully used to re- 
produce the observed column density profiles of dark cloud 
cores (e.g. Alves et al. 120011 Harvey et al. 120011 Harvey et 
al. 120031 Kandori et al. 2005). The BE spheres are described 
in detail by Ebert (119551 1 and Bonnor (I19561 I. As the singu- 
lar isothermal model of S77 the BE spheres satisfy the Lane- 
Emden equation: 



1 d 



fd^ d%> V 

where £ is the dimensionless radius and 
density contrast: 

m = -^-- 

Pc 

The dimensionless radius is linked to the physical radius 
through equation: 
r 

where Ro is the "scale radius" and p c the central density. The 
shape of each BE sphere, i.e. the solution of Eq.[6] can be char- 
acterized by a single parameter, the outer radius of the sphere 
in dimensionless units: % max - The comparison of the observed 
profile with the one derived by solving Eq.[6]thus fixes the pa- 
rameter £ max and the shape of the density distribution. 

Bonnor (1956) examines the stability of the solutions and 
concludes that solutions with £ max value larger than 6.5 are un- 
stable and susceptible for collapse. The solution with £ max = 
6.5 represents the marginally stable configuration. The £ max 
value is directly linked to the density contrast: for example, 
the value £ max =6.5 corresponds to the center-to-edge density 
contrast -14 and the value £ max = 10.0 to ~42. 

While the shape of the density profile is fixed by determin- 
ing the parameter % mdx , matching the scale of the fitted profile 
to the observed profile requires fixing of two additional param- 
eters in Eq. [8] The common way to do this is to fix the sound- 
speed using temperature information derived from e.g. molec- 
ular line observations, and physical radius of the fitted sphere 



{-)(4nG Pc V 



(6) 

is the logarithmic 
(7) 



(8) 
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Fig. 5. Left: The best fitting models of radial density distribution for DC303. The solid line is for a power-law with an exponent 
p = 2.29 + 0.08, the dashed line is a Bonnor-Ebert model with % max = 23 + 3. The dash-dotted line is gaussian convolved average. 
Right: The same for TPN. Solid line is the power-law with an exponent p — 1.19 + 0.09 and dashed lines Bonnor-Ebert models 
with different tj nta x parameter (see text). 



by assuming a certain distance to the globule. However, as dis- 
cussed in Harvey et al. (2003) assuming additional constraints 
such as the outer edge can result in underestimation of the sys- 
tematic errors. A more robust fit is achieved by fitting the scale 
radius, Rq, instead of the dimensionless radius, £ max (see Eq.0. 
In the following we follow this procedure and fit BE models for 
DC303 and TPN. 

The BE model that best fits the shape of the observed col- 
umn density profile of DC303 is described by the parameter 
£, = 23 + 3. As in the case of power-law fitting the inclusion of a 
constant term improves the fit quality significantly. The model 
with a constant (H - K)f g = 0.14 extinction forescreening the 
BE sphere reduces the^ 2 of the fit from 1.60 to 1.47. 

Fitting the scale radius of the BE model did not provide a 
robust fit for TPN, most likely due to non-circular projection 
and the small coverage of the profile compared to the uncer- 
tainties in observed colors. The confusion from the outerparts 
of the globule becomes significant already at ~ 80" distance 
which makes the profile quite noisy. The best fit is marginally 
achieved with a profile having ^ max = 13 + 3, but essentially 
any profile having £ max > 8 produces equally good fit with 
X 1 changing less than one percent. The Fig. [4] shows the fitted 



profiles with % max =8, 12 and 16. Interestingly, the inclusion 
of an additional extinction layer results in a robust fit and, in 
fact, produces similar solution for a density distribution with 
£max = 10 + 2. However, including the forescreen extinction 
does not seem reasonable in this case, as there is no component 
of extinction extending beyond the radius of the BE model. 
Thus we conclude that we can only give lower limit for the 
£max parameter in the case of TPN. 

4.3.3. Masses of the globules 

We calculate the masses of the two globules from the extinc- 
tion maps by summing all the pixels within a certain radius 
together. We adopt the distance of 150 pc derived by Knude & 
H0g (1998) from the data of the Hipparcos satellite. The gas- 
to-dust relation derived by Bohlin et al. d 19781 ) gives the ratio 
of H2 column density and extinction: 



— — = 9.4 x 10 20 cm^mag- 1 . (9) 

Ay 
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Table 1. The summary of fitted models 



Model 




Region 


x 1 


Parameters 




DCld303. 8-14.2 


Single power-law 


24" 


< r < 180" 


1.58 


P = 


2.03 + 0.04 




Single power-law with forescreen 


24" 


< r < 180" 


1.54 


P = 


2.29 ± 0.08, E(H - K) 


= 0.09 


BE 


24" 


< r < 220" 


1.47 




= 23 ± 3, E(H - K) = 


0.14 


Thumbprint Nebula 


Single power-law 


23" 


< r < 100" 


1.13 


p = 


1.19 ±0.09 




BE 


6" 


< r < 120" 


1.11 


£max 


>8±2 





Using the equation above, and adding 30 % due to helium, the 
masses of the cores are: 

M = 5.9 ±1.8 Mo r < 3.4' (0.15 pc) (DC303) (10) 

M = 3.7±1.2M r < 3.4' (0.15 pc) (TPN) (11) 

The masses derived for the cores from the extinction measure- 
ment are in good agreement with the masses derived from the 
C ls O observations (Kainulainen 2004 for DC303; Lehtinen et 
al. [T995l forTPN). 

5. Discussion 

The studies of radial density distributions of dark globules ex- 
ploiting dust emission have established the framework that pre- 
protostellar cores have shallower density profiles and are less 
centrally condensed than protostellar cores. Moreover, the col- 
umn density profiles of the majority of the cores, with or with- 
out central object, seem to be well reproduced with density dis- 
tributions of isothermal Bonnor-Ebert spheres. Despite these 
findings which are made by comparing samples of tens of glob- 
ules, the dust emission is not an ideal tool for characterizing the 
density distributions of individual cores due to interpretational 
problems and limited beam size. 

The extinction measurement using the color excesses of 
background stars is an effort to overcome these uncertainties 
and it provides a direct tracer of mass distribution via high 
signal-to-noise extinction maps. The few previous studies uti- 
lizing the technique in high resolution, this paper included, 
have strengthened the general results of the sub-mm studies: 
the pre-protostellar cores possess density distributions with 
flattened inner regions. 

However, somewhat contradicting the results from sub-mm 
observations, the fitted BE models have often indicated unsta- 
ble configurations even for pre-protostellar cores (e.g. Kandori 
et al. 2005 ). If the starless globules are regarded as stable, hy- 
drostatic configurations, there must be a significant additional 
forces supporting the cloud. Alternatively, the wide range of 
observed BE solutions can be regarded as an evolutionary se- 
quence of slowly collapsing sphere, in which the object is ini- 
tially close to equilibrium and slowly becomes more centrally 
condensed before the actual formation of the central object. In 
this case, the observed central concentration, or the £ max pa- 
rameter, would be directly linked to time since the beginning 
of contraction (Aikawa et al. 2005; see also Fig. 8 of Kandori 
etal. |2005). 



It is necessary to mention that even though the observed 
column density distributions are reproduced using the density 
distribution of the Bonnor-Ebert model, the reality is evidently 
more complex. Even though the extinction features are well re- 
produced, the interpretation of the physical conditions or evo- 
lution of the core is not straightforward nor uniquely deter- 
minable from the extinction alone. The interpretation is ham- 
pered by the unknown contributions from additional supporting 
mechanisms such as magnetic fields and turbulent motions and 
by the basic assumption of thermodynamic equilibrium which 
is not likely met, as well as by deviations from spherical sym- 
metry and non-isothermality. These issues are addressed in de- 
tailed discussion in Hotzel et al. (2002) and Ballesteros-Paredes 
et al. ( 2003). In their work Ballesteros-Paredes et al. argued that 
large fraction of cores which form through turbulent motions 
might resemble the appearance of BE spheres even though they 
are not static but rather dynamically evolving objects. It should 
be emphasized that even though the dust extinction technique 
can trace the column density distribution in high detail, the 
combination of dust emission, extinction and molecular line 
observations is required to make the final arguments about the 
physical properties of the cores. 

5.1. The star forming globule DC303 

In comparison to pre-protostellar globules, which are presum- 
ably close to hydrodynamic equilibrium, it is useful to con- 
sider globules which are evolved to a stage in which the central 
object is established and the system is rapidly accreting mass 
from the envelope (Class YSOs). In the context of dust ex- 
tinction studies, only few occasions which consider cores in 
similar stage can be found. Those are: B335 (Harvey et al. 
1200!) , CB188 (Kandori et al. 12005) and L1014 (Huard et al. 
2006). The column density distributions of these three cores 
can be fitted with highly unstable BE solutions having dimen- 
sionless radii > 15. Particularly, L1014 appears to have a 
very high degree of central concentration with £ mfl x in the range 
29.2 > $ nwx > 36.9. 

The power-law density distribution with p ~ 2 derived for 
DC303 in this study is similar to that of the globule B335, apart 
from the fact that in B335 the evidence of flattening inner re- 
gion at r < 0.03 pc was found. Harvey et al. (120031 ) concluded 
that the observed extinction features of B335 are well repro- 
duced with a broken power-law density structure with infall ra- 
dius of S77 model at 26" ± 3" (0.03 pc at 250 pc), and a simple 
cone-like outflow model. The position of infall radius matches 
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well the molecular line observations and the results of radiative 
transfer modeling following the inside-out collapse framework 
(Zhou et al. 119931 Choi et al. I1995K In these studies the age 
estimate of ~ 1 .3 x 10 5 is derived for the central object. 

If we perform a simple calculation assuming that the shock- 
front of the inside-out collapse is moving outwards from the 
center of DC303 with the speed of a ~ 200 m s _1 , we can 
place an upper limit of ~ 1.8 x 10 5 for the age of the central 
object. Here we assume that we should be able to disentangle 
the flattening from the observations, should it occur at radius 
higher than ~ 50" (0.036 pc). 

B335 is one of the best-known objects harboring a Class 
protostar (see e.g. Hodapp et al. I19981 I. Based on the spec- 
tral energy distribution of the young stellar object in DC303, 
it is presumably in an evolutionary stage where much of the 
mass of the envelope has already accreted to the protostar/-disk 
system. Lehtinen et al. (2005) describes the core to be "a tran- 
sition object between Classes and I". This further indicates 
that DC303 is in a similar stage with B335. The main proper- 
ties of these globules are briefly quoted in table [2] The largest 
difference between the properties seems to be the estimated to- 
tal mass, which is more than double for B335 than for DC303. 

It is interesting that even though DC303 clearly does not 
represent an initial stage of star formation, the normalization of 
the density distribution is well described by the static condition 
of the inside-out collapse. This is in contradiction with studies 
of B335 and L694-2 (Harvey et al. I200T1 120031 where non- 
explainable scaling factors of ~ 3-5 were needed to match the 
scale of the density distribution to the observed color excess. 

The column density profile of DC303 can be equally well 
fitted with a supercritical BE model. The similarity of these 
models originates from the fact that the density of a BE model 
at large radius closely resembles the r~ 2 distribution. At smaller 
radii the profile flattens. As the observed (H - K) colors do not 
show any flattening the BE model produces effectively a power- 
law distribution. Even though the BE model is quite succesfull 
in describing the density distribution, it should be stressed that 
DC303 hardly resembles the configuration of the BE sphere: 
the core has already undergone a collapse forming a central 
object. 

5.2. The starless globule TPN 

Unlike DC303, TPN harbors no embedded source and is in a 
very different stage of evolution. The molecular line observa- 
tions have not shown features indicative for collapse or out- 
flow, but have proven the core to be very cold, T ex « 6 K, and 
close to hydrostatic equilibrium (Lehtinen et al. 119951 [2005 ). 
No significant nonthermal turbulence or large scale motions 
were found in the studies. It should be justified to conclude that 
collapse has not taken place, nor is currently taking place in the 
core. The maximum extinction derived for TPN, A v ~ 13'", is 
also relatively low taking into account that the distribution is 
resolved to a very small radial distance (6", corresponding to 
~ 1000 AU at 150 pc). Correspondingly, the mass of the TPN, 
3.7M , is quite low compared to globules for which the mass 



has been determined using the same method (e.g. Kandori et al. 
120051 Table 5). 

The BE-model of TPN was not conclusively determined, 
but we note the lower limit for the parameter £ mfll > 8. The 
projection of TPN also deviates somewhat from spherical sym- 
metry, which makes the application of BE models more ques- 
tionable. The derived lower limit implies an unstable configu- 
ration disagreeing with the characteristics derived from molec- 
ular line observations. However, this is consistent with the ten- 
dency that a large fraction of globules that are assumed to be 
starless are nevertheless best fitted with unstable BE-models 
(e.g. Kandori et al. 2005). In the context of slowly contracting 
BE-sphere (Aikawa et al. 2005) TPN would be in very early 
stage of contraction with t ~ 5 x 10 5 years since the stable con- 
figuration was actual. There is no detections of infall motions 
in TPN, but the infall speeds at early stages would be too low to 
be detected from the current data. Alternatively, the high cen- 
tral concentration suggests the presence of supporting forces 
in addition to thermal support, and the possibility that TPN is 
not contracting but rather close to equilibrium. With the current 
data it is not possible to conclude whether TPN is slowly con- 
tracting towards the star formation, or close to equilibrium and 
perhaps even dissolving in the future. 

6. Conclusions 

We present deep near-infrared observations of globules DC303 
and TPN made with the ISAAC instrument on the ESO/VLT 
telescope. The H - K colors of the observed stars are used as a 
measure of dust column density and the observed column den- 
sity profile is examined in the light of theoretical models. The 
main conclusions of our work are as follows: 

1. The J, H and Ks band frames of DC303 observed with 
ISAAC show a dense core with peak extinction higher than 
Ay > 50 m surrounding the embedded IRAS source. The 
data of TPN show a more transparent core with peak extinc- 
tion of about Ay ~ 13™. From DC303 the images also re- 
veal Herbig-Haro-like objects close to the core center, and 
a ring-like "halo" structure of scattered light surrounding 
the core. 

2. The best-fit power-law describing the density distribution 
of DC303 has an exponent p = 2.29+0.08 and an additional 
H—K = 0.09 layer of extinction. If the additional extinction 
layer is not included, the power-law exponent has the value 
of p — 2.03 + 0.04. The density distribution shows no signs 
of flattening towards the core center (r > 0.03 pc). The 
normalization of the power-law models is well fitted with 
the static condition of singular isothermal model of S77 and 
the normal reddening law. The density distribution of TPN 
cannot be well fitted with a single power-law. 

3. For DC303, the Bonnor-Ebert model with a dimension- 
less radius £ max = 23 + 3 offers an equally good fit to the 
shape of the observed column density profile as the power- 
law model. The radial density distribution of TPN differs 
greatly from that of DC303 as it shows flattening of the pro- 
file at the distance of about 0.015 pc from the core center. 
We are unable to describe TPN consistently with a Bonnor- 
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Table 2. The properties of DC303 and B335 





DC303 


B335 


Comments 


YSO 


Class 0/1 


Class 


1, 2 


Age 


< 1.8 x 10 5 a 


1.3 x 10 5 a 


3,4 


M 


5.9 MO 


14 MO 


3, from NIR extinction 


"^central object 




0.37 MO 


4 




11 K 


13 K 


5,4 



References: (1) Lehtinen et al. l2005l (2) Hodapp et al. [T998l (3) This 
study (4) Zhou et al. [T993l (5) Kainulainen 120041 



Ebert model, but derive the lower limit for the dimension- 
less outer edge: £ max > 8. 
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